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ABSTRACT 

We report results of an extensivo world-wide observing campaign devoted to a very 
active dwarf nova star - IX Draconis. We investigated photometric behaviour of the 
system to derive its basic outburst properties and understand peculiarities of IX Dra- 
conis as well as other active cataclysmic variables, in particular dwarf novae of the 
ER UMa-type. In order to measure fundamental parameters of the system, we carried 
out analyses of the light curve, O — C diagram, and power spectra. During over two 
months of observations we detected two superoutbursts and several normal outbursts. 
The V magnitude of the star varied in the range 14.6—18.2 mag. Superoutbursts oceur 
regularly with the supercycle length (P sc ) of 58.5±0.5 d. When analysing data over the 
past 20 years, we found that P sc is increasing at a rate of P = 1.8 x 10~ 3 . Normal out- 
bursts appear to be irregular, with typical oceurrence times in the range 3.1 — 4.1 d. 
We detected a double-peaked structure of superhumps during superoutburst, with 
the secondary maximum becoming dominant near the end of the superoutburst. The 
mean superhump period observed during superoutbursts is P s h = 0.066982(36) d 
(96.45 ± 0.05 min), which is constant over the last two decades of observations. Based 
on the power spectrum analysis, the evaluation of the orbital period was problem- 
atic. We found two possible values: the first one, 0.06641(3) d (95.63 ± 0.04 min), 
which is in agreement with previous studies and our O — C analysis (0.06646(2) d, 
95.70 ± 0.03 min), and the second one, 0.06482(3) d (93.34 ± 0.04 min), which is less 
likely. The evolutionary status of the object depends dramatically on the choice be- 
tween these two values. A spectroscopic determination of the orbital period is needed. 
We updated available information on ER UMa-type stars and present a new set of their 
basic statistics. Thereby, we provide evidence that this class of stars is not uniform. 

Key words: binaries: close - stars: cataclysmic variables, dwarf novae, individual: 
IX Draconis 



1 INTRODUCTION 

Cataclysmic variable (CV) stars are close binary systems 
which consist of two component stars: a donor (secondary) 
star which fills its Roche lobe and loses its mass through the 
inner Lagrangian point onto a white dwarf (primary star). 
Dwarf novae are one of the non-magnetic types of CV sys- 
tems. They have a white dwarf as a primary, and typically 
a low-mass main sequence star as secondary. The magnetic 
field of the white dwarf is weak enough (< IO 6 G), so that 
the material from the secondary is transported onto the pri- 
mary through an aceretion disk around the white dwarf. 
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Such a disk is the m ain source of optical radiation of the 
system l|Hellierll2003 h Mass transfer in dwarf novae leads to 
outbursts which are explained by the thermal-tidal instabil- 
ity model (TTI) as caused by a sudden gravitational energy 
release due to aceretion of material onto the white dwarf 
(|Osakilll974l h 

One of the types of dwarf novae are SU UMa stars which 
are characterized by short orbital periods (P or b < 2.5 h). 
Their most important feature is that, beside normal out- 
bursts, they show superoutbursts about 1 mag brighter and 
about ten times less frequent than normal outbursts. During 
superoutbursts there are observed so-called superhumps, i.e. 
additional "tooth-shape" modulations in the light curve with 
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amplitudes of a fraction of a magnitude, and periods larger 
by about 3% than the corresponding orbital period. 

The most active a mong them are ER UMa-type stars 
l|Kato fc Kuniavalll995h . These are objects with extremely 
short supercycle lengths of about 20 — 60 d, normal outburst 
cycle lengths of about 3 — 4 d and amplitudes of outbursts 
of about three magnitudes, which is small in comparison to 
typical SU UMa stars. They seem to be standard SU UMa 
stars with only a larger act ivity and lum inosity due to their 
higher mass transfer rates (|Osakilll996l ). 

The length of the supercycle and the number of normal 
outbursts in between superoutbursts are specific for each 
SU UM a-type star. This feature can be explained by the TTI 
model (Osaki 1996). According to this model, normal out- 
bursts remove less material from the disk than was accreted 
since the last outbursts. Thus, the disk is growing during 
successive normal outbursts. After it reaches the criticai ra- 
dius (r ~ 0.46a, with a being the binary separation), the 
disk becomes eccentric and tidal effects occur. This leads to 
a superoutburst, during which the disk is shrinking again. In 
this scenario the sypercycle length is set by the mass transfer 
rate in the way that the observed length of the supercycle 
is inversely proportional to the mass transfer rate. Accord- 
ing to the standard theory of evolution of CVs, an object 
evolves from long to short orbital periods and from high 
(~ 1 0~ 8 M B /yr) to low (~ 10" 10 M /yr) mass transfer 
rates (|Knigge et al.ll201ll ). That is why it is very surpris- 
ing to find extremely active ER UMa-type stars among the 
objects with the shortest orbital periods, or even as period 
bouncers (which should have the lowest mass transfer rate 
~ IO" 12 M /yr). 

Until now there were discovered only five stars be- 
longing to the ER UMa-type of dwarf novae, among them 
IX Draconis (IX Dra). Thus, studying each one of them is 
of primary importance. 

IX Dra was detected for the first time bv lNoguchi et ali 
(1980), and described as a variable object wit h an UV excess 
and a proper motion of 0.016 arcsec/yr by iNoguchi et al.l 
(|1982T) . The o bject was ide ntified as a cataclysmic variable 
candidate by iKlosd (I1995Í ). based on the photometric be- 
haviour and the fact that it has an UV excess. IX Dra 
was mentione d in the "Cataloq a nd Atlas of Cataclysmic 
Variables" by iDownes et ali fl997) as a dwarf nova candi- 
date of the U G em class, variabl e in the magnitude range 
of 15.0 — 18.5 B. iLiu et al.l |l999) published a spectrum of 
IX Dra obtained at V magnitude of 16.3 which showed a 
strong blue continuum and emission lines indicating that 
the spectrum was not taken during quiescence. This was the 
definite proof of the CV nature of this objec t . A d etailed 
study on IX Dra was reported bv llshioka et al.l (|200lh . They 
performed a few-months-long observational campaign and 
found that the object belongs to the ER UMa-type of stars. 
They were able to estimate the basic observational parame- 
ters of the system: supercycle length of 53 d, normal outburst 
cycle length of 3 — 4 d, outburst amplitude of 2.5 mag, 
and suspected superhump period of 0.067 d. Another ex- 
te nsive photom e tric o bservational campaign was presented 
bv lOlech et al. | <|2004l ). They reported similar estimates of 
the main observational parameters of IX Dra: supercycle 
and cycle lengths equal to 54 ± 1 d and 3.1 ± 0.1 d, respec- 
tively. Additionally, they estimated the superhump period as 
0.066968(17) d and found another period of 0.06646(6) d in 



the light curves. If the latter is the orbital period, the value 
of the superhump period excess would suggest a very low 
mass ratio of the binary system, which in turn implies that 
IX Dra has a brown dwarf as secondary. Tha t would make 
IX Dr a one of the most evolved dwarf novae. Ilshioka et alj 
(2007) presented results of infrared spectroscopy which were 
performed in order to determine the spectral types of sec- 
ondary stars of a few CV systems including IX Dra. For 
this system, they found typical features of late-type dwarfs, 
with the best fit to the spectrum corresponding to a tem- 
plate spectrum of M1-M5 dwarfs with a power-law compo- 
nent with an index —3.1. What is more, from the relation 
between the absolute magnitude and the period, they esti- 
mated that the distance to the object is 995 pc, and from 
the template dwarfs of M1-M 6 types the expected distance 
is 1000 - 3000 pc. Although lfshioka et alj <|2007h consider 
that the secondary is a more massive star, the possibility 
that the secondary is actually a brown dw arf with a large 
mass- transfer rate could not be ruled out. iPretorius et al.l 
(2007) estimated the distance as d = 430ÍfgQ pc based 
on the apparent brightness in the infrared. These recent 
facts were our motivation to study this interesting object 
and to verify its evolutionary status. It is very important in 
the context of our understanding of the still poorly studied 
ER UMa-type of stars. 

The paper is arranged as follows. In Section [2] we give 
information on our observations and data reduction. The 
global light curve of IX Dra is presented and analysed in Sec- 
tion [3] In Section [4] we analyse the periodicity of detected 
superhumps. Section [S] is dedicated to a wider interpreta- 
tion of properties of IX Dra and the rest of members of the 
ER UMa-type class. We summarize the conclusions of this 
work in Section [5] 



2 OBSERVATIONS AND DATA REDUCTION 

Observations of IX Dra reported here were obtained between 
2010 October 3 and December 13. During this time the dwarf 
nova showed two superoutbursts (the first one was observed 
during 15 nights, and the second one during 11 nights) and 
seven normal outbursts (observed during 34 nights) . The ob- 
servations in total cover 60 nights with 10613 measurements 
with the average exposure time of about 120 s, depending 
on the instrument. 

Altogether IX Dra was monitored by six observers us- 
ing eleven telescopes with primary mirror diameters rang- 
ing from 0.25 to 1.0 m, located in the USA and spread 
over Europe. Such locations of telescopes allowed us to 
obtain a beautiful light curve without significant breaks, 
especially during superoutbursts which are of the great- 
est importance for us. For much of that time observa- 
tions were performed almost without breaks, the main lim- 
itation came from the weat her. The observer s are mem- 
bers of the CURVE team (jOlech et alj |200S| . 120091 ). see 
http : //users . camk.edu.pl/maRdaot/ curve/, and the Cen- 
ter fo r Backyard Astrophysícs ( CBA ) (|Patterson et al.ll2003Í . 
2005), see http://cbastro.org/. 

The first of our instruments used during this observa- 
tional campaign was the 0.6 m telescope of the Warsaw Uni- 
versity, located in Ostrowik, Poland. The observations were 
also carried out at the TUBITAK National Observatory in 
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Table 1. Our obscrvations of IX Dra in 2010. 



Observatory 


Country 


Telescope 


Observer 


Warsaw University Observatory 


Poland 


0.6 m 


M. Otulakowska-Hypka, A. Olech 


TUBITAK Observatory 


Turkcy 


0.4 m, 0.6 m, 1.0 m 


A. Rutkowski 


Observatório dei CIECEM 


Spain 


0.25 m 


E. de Miguel 


Antclope Hills Observatory 


USA (CO) 


10" 


R. Koff 


NF / Observatory 


USA (NM, AZ) 


0.6 m, 24" 


A. W. Neely 


Tzec Maun Foundation Observatory 


USA (NM) 


14", 16", 18 cm 


K. Bakowska, M. Otulakowska-Hypka 



Turkey with 0.4 m, 0.6 m, and 1.0 m telescopes. IX Dra was 
also observed remotely using small telescopes of the Tzec 
Maun Foundation, located in New México, USA. Additional 
observations were done with the Webscope (24") and Su- 
per Lotis (0.6 m) telescopes of the NF/ Observa tory in New 
México and Arizona, USA (jCrauer et al.ll200Sl ). The input 
from the CBA carne from a 0.25 m telescope at Observatório 
dei CIECEM, Huelva, Spain, and a 10" telescope at Ante- 
lope Hills Observatory, Colorado, USA. Information on the 
observations are gathered in Tab. [T] 

The object was observed with the clear filter ("white 
light") in order to obtain precise photometry of the star in 
quiescence and also at minimurn light - below 18 mag. Flat- 
field, dark, and bias corrections were performed on the raw 
files in a standard way. In the case of the observations with 
small telescopes from CBA, this was done with commercially 
available software. For the rest of the observations, standard 
IRAF routinesE] were used. To extract the pro file photometr y 
we used DAOphot II and Alistar packages |Stetsonll 19871 ). 
Relative unfiltered magnitudes of IX Dra were derived from 
the difference between the magnitude of the object and the 
mean magnitude of three comparison stars which are shown 
in Fig. □ 

Instrumental magnitudes were transformed into 
the standard V sy st em u sing photometric data from 
lHenden fc Honevcuttl |l995). Such a transformed global 
light curve of IX Dra spanning over two months of 
observations is presented and analysed in the next section. 

Typical photometry errors of our measurements are in 
the range 0.007 — 0.21 mag during the bright state, and 
0.01 — 0.25 mag at the minimurn light. The median value of 
the precision equals to 0.019 and 0.023 mag, respectively. 



3 GLOBAL LIGHT CURVE 

For convenience, from now on we use only a day number, 
HJD — 2455000 [d], to refer to our observations. Our observ- 
ing campaign started on day 473. 

Figure [2] presents the light curve of IX Dra during the 
whole campaign. We were lucky to start the observations at 
the time of appearance of the precursor of the first of two 
superoutbursts. On day 476 the superoutburst started and 
lasted for 15 days (2 days of initial rise, 11 days of the plateau 
phase, and 2 days of final decline). The coverage of this part 

1 IRAF is distributed by the National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under a cooperative agreement with 
the National Science Foundation. 
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Figure 1. Position of IX Dra and three comparison stars on a 
chart. 



of the light curve turned out to be the best of the campaign, 
as we were able to collect data during every single night from 
at least one of the observatories. This superoutburst was 
followed by a set of normal outburst which in total lasted 
for 43 days. We were only able to observe IX Dra on 31 of 
these nights, due to adverse weather conditions. On day 534 
the second superoutburst started and was observed over the 
following 10 days. Weather conditions were not as good as 
during the first superoutburst. 

The most immediate conclusion from the light curve is 
probably the fact that IX Dra is an extremely active dwarf 
nova. Every outburst is almost at once followed by another 
one. There is no evident quiescence period visible in the light 
curve between these outbursts. This photometric beh aviour 
is typ ical for dwarf novae stars of the ER UMa-type (|Osakil 
119961 ). 

During superoutbursts the brightness reaches V ~ 
14.7 mag at maximum, and fades to V «18 mag 
at minimurn, resulting in the superoutburst amplitude of 
A s ~ 3.3 mag. In turn during normal outbursts the bright- 
ness range spreads from V w 15.3 mag at maximum 
to V ~ 18.1 mag at minimurn, giving the amplitude of 
A n ~ 2.8 mag. These values are typical for ER UMa stars, 
see Tab. H 
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Figure 2. Global light curve of IX Dra during our observational campaign 
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3.1 Supercycle length 

Until now there are only three previously reported mea- 
surements of the supercycle length (Ps c) for I X Dra . The 
first one, P} c = 45.7 d, was given by iKlosd |l995) with 
no reported estimate of the corresponding uncertainty. The 
author had very unevenly sampled data, so the periodicity 
analysis was challenging. This uncertainty, however, can be 
significant, as the author claimed that the period P 3C did not 
match ali his dat a. Also the second m easurment, P 3C = 53 d, 
was presented bv llshioka et alJ l)200lf ) with no uncertainties. 
In this case authors did not give a reason for that. The third 
measurement of the supercycle length, P 3 C = 54 ± 1 d, was 
presented bv lOlech et al.l |2004h . Considering the probably 
large uncerta i nty o f the first of these value s, the results of 
llshioka et ãH (120011 ) and lOlech et al.l l|2004h seemed to sug- 
gest a constant value of P sc for IX Dra. 

To derive P sc from our observations, we com- 
puted the power spe c trum using the ANOVA code of 
ISchwarzenberg-Czernvl (|l996T l only for the part of the light 
curve which covers both superoutbursts. We found that the 
most prominent peak of the power spectrum is the one at 
the frequency 0.0171(2) c/d i.e. P sc = 58.5 ± 0.5 d. 

Thus, we provide a new estimate of the supercycle 
length for this object and thereby we show that P 3C has 
increased during the past twenty years (see Fig. [3j . The cor- 
responding rate of increase of the period is P — 1.8 x IO -3 . 

In spite of the fact that our estimate of P 3C is based 
on only two observed superoutbursts, such a change of the 
supercycle length seems reasonable, because such changes 
were also observed for other ER UMa-type stars. For in- 
stance lKatd (|200lh reported an extremely variable supercy- 
cle length of V1159 Ori. In his Fig. 2 one can see rapidly 
changing jumps of this value between 44.6 and 53.3 d over a 
span of about five years. RZ LMi has a supercycle length 
decreasing with the rate of P = —1.7 x IO -3 , in turn 
ER UMa has t he supercycle length in creasing with the rate 
P~4x IO" 3 (|Robertson et al.lll995l l. 
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Figure 3. The increasing supercycle length of IX Dra during the 
past twenty years. The data points are taken from the literature. 
Uncertainties are given when available. The line represents the 
best fit to the data. The corresponding rate of increase of the 
period is P = 1.8 X IO -3 . 



We used the P sc to create the folded light curve of our 
both superoutbursts. It is shown in Fig. [4] 

Although the light curve folded with the P 3C looks rea- 
sonable, and despite the fact that significant changes of P sc 
are apparently typical for ER UMa-type stars, we decided to 
investigate our result more precisely to make sure that it is 
not an incidental effect based on only two perhaps untypical 
superoutbursts. We looked for additional recent photomet- 
ric data of the object, which cover a greater time range. In 
the archive of The American Association of Variable Star 
Observers (AAVSO) we found a light curve for IX Dra. Its 
time coverage was quite poor, probably because of the fact 
that the object reaches 18 mag in quiescence and thus it 
is then too faint for amateur observers. However, we were 
lucky to find one superoutburst with sufncient time cover- 
age which was just following our observational run. In Fig. [5] 
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Figure 5. Observed superouturst of IX Dra (points) from our observational campaign and the succeeding superoutburst from the 
AAVSO archive, suggesting that the P sc = 58.5 d is consistent with ali three superoutbursts. Lines indicate anticipated occurrences of 
superoutbursts with the supercycle length P sc . 
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Figure 4. Light curve of IX Dra during superoutbursts obtained 
by folding the light curve with P BC = 58.5 ± 0.5 d. Grey circles 
indicate the data of the first superoutburst, black diamonds rep- 
resent the data of the second superoutburst. 
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Figure 6. Power spectrum for the light curve of normal outbursts 
with the highest peak corresponding to P~ = 4.1 ± O.ld. 



a brown dwarf donor, a nd that the high m ass transfer rate 
is an unusual state (see lOlech et al.l (|2004h and Sec. 14.31 for 
details). 



we present both superoutbursts observed during our cam- 
paign together with the succeeding superoutbursts from the 
AAVSO archive (at day ~ 600) which is in agreement with 
our supercycle estimate P sc ~ 58.5 ± 0.5 d. The superout- 
burst profile, which is shown as the solid line in the Fig. O 
was obtained by completing an analytical fit using Bezier 
curves to the first superoutburst. The derived profile was 
copied forward with supercycle period P 3C . 

There is one more interesting issue possible to check 
out in the context of the increase of the supercycle length of 
IX Dra: the future evolution of the system. Assuming that 
the supercycle length will increase in the same way in the 
future, we can predict when IX Dra would reach the status 
of SU UMa and WZ Sge stars. Typical supercycle lengths for 
these types of stars are a few hundred days and decades, re- 
spectively (|Hellierll200Íh . Thus, assuming that this increas- 
ing trend of P sc will continue, we calculated that IX Dra 
would attain the stage of SU UMa with the P sc = 300 d in 
almost 400 years from now. In turn in about 7000 years the 
system would become a WZ Sge-type star with the super- 
cycle length of 10 years. We wish to stress that this is only 
a speculation which is in contradiction with the most likely 
scenario of the evolution of the object i.e. that IX Dra has 



3.2 Normal cycle length 

In contrast to the supercycle, there is not such a straight 
periodicity in the observed occurrences of normal out- 
bursts. We computed power spectra for the light curve cov- 
ering only normal outbursts , using the ANOVA code of 
ISchwarzenberg-Czernvl |l996). We ran a number of tests 
with a variety of frequency ranges and steps but without any 
fully satisfactory result. One can say that the most promi- 
nent peak is the one shown in Fig. [5] and corresponding to 
the frequency 0.24 ± 0.05 c/d i.e. P C J = 4.1 ± 0.1 d. 

However, we cannot accept this value as the normal cy- 
cle length of the system, P c . The folded light curve with such 
a period did not give an adequate result. Our conclusion is 
that normal outbursts observed during our campaign did 
not occur at regular intervals. On the other hand, we know 
the previousl y estimated value of the normal cycle length of 
IX Dra from íOlech et al.l (I2004T ): P c 2 = 3.1 ± 0.1 d. To make 
this issue clear, we plot the light curve with superimposed 
two curves indicating regular anticipated occurrences of nor- 
mal outbursts with Pç and another one with P*. First, we 
made a fit to the last ouburst (at 531 d), which is the only 
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HJD-2455000 [d] 

Figure 7. Normal outbursts of IX Dra during our observational campaign do not appear to be periodic. We show the observational 
data (black dots) with two superimposed curves indicating regular anticipated oceurrences of normal outbursts with the most likely 
period estimated from our observations Pç = 4.1 d (black dashed line), and with the value of the normal cycle length from the literature, 
P£ = 3.1 d (grey solid line). Neither period decribes the lighteurve well. 



one with a sufiicient time coverage, and then repeated both 
fits backward in time with the given periods. See Fig. [7] 

It is now evident that normal outbursts are not exactly 
periodic, their oceurrence seems to be pretty random. Only 
two of them, observed at days ~ 506 and 531, emerged more 
or less according to both predictions, Pç and Pç. 



4 SUPERHUMPS 

The light curve of IX Dra shows clear superhumps during 
superoutbursts as well as during normal outbursts. In Fig. [8] 
we show nightly curves of the first superoutburst to illustrate 
their behaviour. 

The first three nights (473 — 475) correspond to the 
precursor and show small amplitude modulations of about 
0.15 mag. The next two nights (476 — 477) correspond to the 
quick initial phase of the superoutbursts and show modula- 
tions with the amplitude reaching 0.2 mag. Just after hitting 
the maximum brightness of the light curve (days 478 — 479) , 
there are visible clear and aceurate, so-called "tooth-shape", 
superhumps with the maximum amplitude 0.35 mag, de- 
creasing with time to 0.24 mag during the whole plateau 
phase. A very interesting behaviour begins at day 480: aside 
from the normal superhumps we observe secondary máxima 
shifted by 0.5 in phase with respect to primary superhumps, 
with the amplitude of about 0.05 mag increasing with time. 
After day 482 the light curve gets complex and even chaotic 
- it is quite tricky to tell apart the primary and secondary 
superhumps. During this time, the amplitude of primary su- 
perhumps is decreasing and the amplitude of the secondary 
humps is increasing. At the end of the superoutburst (day 
490) we see again strong clear humps, with the amplitude 
reaching 0.4 mag. Although, their shapes differ from the 
typical "tooth-shape" of classical superhumps. 

To investigate the apparent periodicities in the light 
curve, we performed power spectra and O — C analyses. They 
are described in the following subsections. 



4.1 Period analysis 

At first, to refine our period search, we converted the light 
curve in two ways. For a start, we transformed the mag- 
nitudes into intensities to equalize the amplitudes of the 
superhumps. Then, we detrended the light curve for each 
day separately, by subtracting a best-fitting parábola. As a 
result we obtained light curves in units of counts with the 
average intensity at zero. 

Due to the differing quality of the data of each super- 
outburst, we decided to analyse them separately. We divided 
the whole data into three data sets, for: 

(i) the lst superoutburst, 

(ii) the 2nd superoutburst, and 

(iii) normal outbursts and quiescence. 

For each of these subsets, we performed the period anal - 
ysis with the ANOVA code of lSchwarzenberg-Czernvl (|l996l ) 
to calculate power spectra. Results are presented in Fig. 

4-1.1 Superoutbursts 

For both superoutbursts there is clearly visible one promi- 
nent peak corresponding to the superhump period, P s h- 
For the first superoutburst it is located at the fre- 
quency 14.910(7) c/d which determines the period P} h = 
0.06707(3) d, and for the second one the main peak lies at 
the frequency 14.96(2) c/d, i.e. P^ h = 0.06683(9) d. These 
period uncertainties are 1 sigma uncertainties. 

We adopt P} h as the subsequent superhump period. It is 
more reliable than P^ h due to a greater time coverage of ob- 
servations of the first superoutburst and consequently much 
better S/N ratio of its periodogram. We can say that the 
result from the second superoutburst, taking into account 
its worse quality, agrees with P} h to within 2 sigma. 

We now consider the analysis of the stability of the 
phase of superhumps. The question about the presence of a 
phase reversal in ER U Ma-type star s was t aken up by a num- 
ber of researchers, e.g. iKato et al i (|2003l ) and lOlech et al.l 
(2004). To check the stability of the phase of superhumps, 
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Figure 8. Supcrhumps at each day of the first superoutburst. Consecutive nights are denoted by dates given as HJD- 2455000 [d] 



8 M. Otulakowska-Hypka, A. Olech et ai 



600 
500 
400 
300 
200 
100 




180 
160 
140 
120 
100 
80 
60 
40 
20 














































-li 














- — ■ ' -1 


K.J..I, .. . 



10 



15 



20 25 
Frequency [c/d] 



30 



35 



10 



15 



20 



25 



30 



35 



Frequency [c/d] 




10 15 20 25 30 35 

Frequency [c/d] 



Figure 9. Power spectra for the light curve (from top to bottom) 
of the lst superoutburst, the 2nd superoutburst, normal outbursts 
and quicsccncc. 



at first we removed the general decreasing trend from the 
light curves, and then phased them with the superhump pe- 
riod for each of the nights separately. In Fig.[TD]there is such 
a plot for the data of the first superoutburst. For each of the 
light curves we also plot a corresponding fit to the phased 
data. 

It is now clearly visible that during the first two nights 
we observed only the primary superhumps with máxima at 
phase ~ 1. On day 480 we detected additional humps in the 
light curve which begins to appear at the phase ~ 0.5. It 
is present as from now, and its amplitude is continuously 
increasing over the next days. On day 485 we can see that 
the amplitude of the secondary humps gets higher than the 
amplitude of the primary superhumps, which is gradually 
decreasing over the time. The phase of the secondary humps 
is also changing - it is decreasing from ~ 0.5 to ~ 0.1. 
Finally, at the day 490, the place of the primary superhumps 
is taken by the secondary humps. As was mentioned before, 
the shape of these final humps differ from the initial "tooth- 
shape" form. 

We performed prewhitening of the detrended light curve 
in order to search for additional periodicities. At first, from 
the light curve we removed a modulation with the super- 
hump period, which is shown in Fig. [T0]as a black curve. On 
such data we performed the power spectrum analysis. The 
result is presented in Fig. 1111 

The power spectrum is pretty noisy, although one could 
claim that there is some signal around the frequency ~ 9 c/d, 
and probably its first harmonic at ~ 20 c/d. However, this 
would be unexpected and unlikely. We detected no trace of 
orbital modulation during the superoutburst. As reported 
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Phase 

Figure 10. Evolution of superhumps during the first superout- 
burst. Plots show the detrended light curve (grey points), phased 
with the superhumps period, P\, for each day separately. Num- 
bers on the right indicate dates of observations (HJD-2455000 
[d]). The black curve represents our fit to the light curve. See the 
tcxt for more details. 
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Figure 11. Power spectrum for the light curve of the first super- 
outburst prewhitened with P B h- 
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Figure 12. Zoom in the power spectrum of the light curve with 
normal outbursts showing a set of main peaks around / = 15 c/d. 



frequencies (c/d) periods (d) 



h 


= 13.990(7) 


Pfi 


= 0.07148(4) 
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= 14.826(7) 




= 0.06745(3) 


h 


= 15.057(7) 




= 0.06641(3) 


h 


= 15.428(7) 


Pf, 


= 0.06482(3) 



Table 2. The most prominent frequencies for the periodogram of 
normal outbursts. 

by lOlech et ali (|2004Í ). the orbital signal for IX Dra was 
observed duringthe superoutburst in September 2003. Thus, 
we also expected to detect it at the frequency ~ 15.05 c/d. 

4-1.2 Normal outbursts and quiescence 

When it comes to the power spectrum of the light curve 
covering normal outbursts and quiescence, it does not show 
definhe evidence of periodicities (see again Fig. O. However, 
one can distinguish a number of peaks with similar power 
around ~ 15 c/d, corresponding to Pj h and P^ h . These most 
prominent frequencies are labeled as /i, /i, /3, and fi in 
Fig. 1121 where we show the power spectrum zoomed into 
the region of interest. The associated periods are presented 
in Tab.[2] An interpretation of these signals will be deferred 
to the following sections. 

4.2 The O - C analysis 

IX Dra shows clear humps during entire superoutbursts and 
also throughout normal outbursts and the very short qui- 
escent state. In the previous subsection we found a clear 
periodicity in the data of both superoutbursts, and an indi- 
cation of the presence of a similar periodicity between them. 

Previous studies of this object indicate that apart from 
typical superhum ps, there is an add itional light curve mod- 
ulation observed joiech et al.ll2004T ). 

To perform the O — C analysis of our data, we care- 
fully determined moments of ali máxima and, as before, we 
examined them for each part of the light curve separately. 

4-2.1 Superoutbursts 

We present results based on only the first superoutburst. 
This data set is more credible because of its much better 
quality and full time coverage, in contrast to the second 
superoutburst. 

In this part of the light curve we identified 60 máxima 
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Figure 13. Light curve of the first superoutburst of IX Dra (top 
panei), evolution of the amplitude of superhumps (middle panei), 
and O — C diagram for corresponding superhumps máxima. In 
the two lower plots black dots indicate typical integer values of 
cycle numbers, open circles present values shifted by a half of a 
cycle. Details in the text. 



of superhumps. A linear fit to the data brought the following 
ephemeris: 

HJD max = 2455478.2865(5) + 0.066894(6) x E (1) 

The corresponding O — C diagram is presented in 
Fig. El 

We list times, errors, cycle numbers, and amplitudes for 
each superhump in Tab.(3] There are also given the O—C val- 
ues calculated according to the obtained ephemeris (Eq. [T]). 

There are two very intriguing things about this O — C 
investigation. The first peculiar issue is the fact, that be- 
side standard integer values of cycle numbers, we found 26 
máxima which are shifted by a half of a cycle with respect 
to these typical ones. They are shown in Fig. [13] as open 
circles. Another interesting matter is the jump of the phase 
by about 0.15, observed somewhere around the cycle num- 
ber E — 60. Before and after the jump ali the data points 
arrange as inclined straight lines. 

We separated the data into two data sets (E < 60 
and E > 60) in order to check out if we observe only 
the phase jump itself, or rather a period change. We fitted 
linear functions to both data sets and obtained the following 
ephemerides for E < 60 and E > 60, respectively: 

HJD max = 2455478.2865(7) + 0.066841(23) x E (2) 
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Table 3. Cycle numbers (E), times of máxima (T max , i.e. HJD — 
2455000 [d]), O — C values, and amplitudes (A) of supcrhumps 
observed during the first superoutburst. 
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HJD max = 2455478.306(2) + 0.066746(14) x E (3) 

The difference between these two periods is equal to 
0.000095 d, which is four times greater than the accuracy 
of measurements. This result is quite ambiguous because it 
is almost bound together with the 3 sigma extent of the 
periods. Thus, a conclusion that at about the cycle E ~ 60 
there exists a change between two separate periods would 
probably be exaggerated. 

We claim that this is the effect of only a phase shift 
which is also often observed in other SU UMa stars. An 
explanation of this feature is that at some point the main 
source of the super hump signal starts to be dominated by th e 
bright spot rather than the fiexing disc (|Wood et al.ll201lT ). 
This is expected to happen at the moment of transition from 
the plateau phase to the decline of a light curve. However, 
in the case of this particular superoutburst of IX Dra, the 
phase shift occurs before the transition, which is unusual. It 
is unclear if this feature is specific to this particular super- 
outburst or it is an inherent characteristic of the system. Pre- 
viously published light curves of IX Dra, which show double 
periodicity, are too complicated to identify secondary humps 
|Qlech et al.ll2004D . 

On the other hand, in another ER UMa-type object, 
DI UMa, secondary máxim a are detected only on the very 
last days of superoutburst (fRutkowsk i et aÍ1l2009l . Fig. 7). 
It is also worth mentioning that the occurrence of such sec- 
ondary máxima, shifted by a half of a cycle with respect 
to primary superhumps, is most probably a characteristic of 
on ly dwarf novae with shor t orbital periods. It was suggested 
bv ÍRutkowski et al.l |2007Í ) that objects with longer orbital 
periods are less likely to show such behaviour. 

4-2.2 Normal outbursts and quiescence 

We performed the O — C analysis also for the part of the 
light curve between both superoutbursts, covering a set of 
normal outbursts and a quiescent state. Although the time 
coverage here was much poorer, we were able to detect a 
number of clear humps. 

Modulations of the light curve around máxima of nor- 
mal outbursts are very weak. There is only a subtle flickering 
there, like in the case of typical SU UMa stars. Clear humps 
were identified only below the mean magnitude of outbursts, 
mainly very close to the minimum brightness of the star. 

The power spectrum for this part of the light curve, 
shown in Sec. 14.1.21 is very complex in structure and shows 
several parallel peaks around the value / ~ 15 c/d. This 
suggests only the fact that there is some periodicity present, 
but without any quantitative estimate. That is why it was 
crucial to perform as optimal an O — C analysis as possible, 
to verify the accurate value of a period observed during nor- 
mal outbursts and quiescence. Thus, we precisely measured 
times of máxima of every single hump which we were able 
to catch, even those with a lower amplitude, or with a poor 
time coverage of observations. 

The O — C diagram for ali these points is rather noisy. 
Nonetheless, a linear fit to these points allowed us to deter- 
mine the ephemeris: 

HJD max = 2455497.73(1) + 0.06639(3) x E (4) 



An interesting thing is the fact that the obtained period 
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Table 5. Pcriods derived from our data. Values are expressed in 
days. 
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Figure 14. The O — C diagram for máxima of definite humps 
observed during normal outbursts and the quiescent state. 



Table 4. Cycle numbers E, times of máxima (T ma x, ic HJD — 
2455000 [d]), and O — C values (in units of days) of clear and def- 
inite humps observed during normal outbursts and the quiescent 
state. 
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0.06639(3) d is in agreement with the P/ 3 , which in our 
opinion could be the orbital p eriod of the s ystem , and with 
the orbital period identified in lOlech et al. (2004). However, 
this data set does not indicate any information about the 
stability of the obtained period. 

To check the stability of this period, we decided to per- 
form one more improvement and restrict our O — C analysis. 
From the previous data set of ali possible humps, we sepa- 
rated only 6 humps which are obviously clear and accurately 
covered by observations, thus are beyond any doubt. For this 
sparse data set we performed another O — C analysis. 

From a linear fit to these points we obtained the follow- 
ing ephemeris: 

HJD max = 2455508.176(3) + 0.06646(2) x E (5) 

The consequent O — C diagram is shown in Fig. 1141 
Times, errors, number of cycles (E), and O — C values 
computed according to the linear ephemeris for each of these 
points are given in Tab. [4] 

Although based on this very limited data set, we ob- 
tained a result which is in agreement with the orbital period 
Pf 3 (Eq. [5]l . What is more, the stability of this value is now 
confirmed on the basis of the O — C diagram (Fig. I14|l . 



4.3 Periods' interpretation 

AU relevant periods emerging from our photometric data 
have been collected in Tab. [5] 

Periods obtained during superoutbursts are consistent 
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between both methods of analysis. We detected only one 
clear periodicity, which is the superhump period, P s h- As 
a representative value, we consider here the corresponding 
mean value, namely, P s h — 0.06698(4) d. 

The obtained P 3 h indicates that the value of the super- 
hump peri od remains constant for t his object over the last 
decade, see llshioka et al.l (|200ll ) and lOlech et al.l (|2004| ). Al- 
though, in contrast to the latter publication, we did not 
detect any orbital period modulation of the light curve 
throughout superoutbursts, which is surprising. 

The presence of an additional hump during superoutu- 
bust (see Fig. [TÕ] and Fig. 1131) has been previously dis- 
cussed, but it requires interpretation. It appears just after 
four days since the beginning of the superoutburst (about 
the cycle number E ~ 60) with the same period as typi- 
cal superhumps (P s h) but with a half-cycle phase shift. At 
about E ~ 135, amplitudes of these secondary humps be- 
come dominant relative to primary superhumps. As a conse- 
quence, we observe a jump of the phase in the O — C diagram 

(Fig.na. 

An explanation of this feature could perhaps be the 
presence of an additional source of flux on the surface of 
the accretion disk of IX Dra. It is worth mentioning here 
an altern ative interpre tation of superhumps which was pre- 
sented bv ISmakl (|2007l ) when analysing a SU UMa-type ob- 
ject, Z Cha. He showed that it is possible that the stream 
of matter overflows the disk of a system, and beside the 
"standard" hot spot, it causes additionally a second "mod- 
ified" spot, presumably on the other side of the disk. Per- 
haps this could be the source of such double humps? Such 
an explanation is tempting and it would imply a need 
of an improvement of the standard, c ommonly accepte d 
tidal-resonanc e mod el of superhumps |Whitehurstl 1988). 
lUemura et al.l (j2012) also showed a model with two flaring 
parts in the disk of SU UMa-type star, which produce early 
superhumps with secondary humps. They also claim that the 
tidal effects alone canno t expl ain such a behaviour. On the 
other hand, IWood et al.l (|201ll . and references therein) men- 
tioned an example of a "double humped" light curve with 
two approximately antiphased signals, and show a model 
with two distinct physical mechanisms which could cause 
positive superhumps, namely the viscous dissipation in the 
flexing disk and the time-variable viscous dissipation of the 
bright spot. Also, a proposal of a new interpretation of su- 
perhumps as they are caused due to enhanced dissi pation 
of the kinetic energy of the stream was suggested bv ISmakl 
2009). This effect is a consequence of a modulated irradi- 



ation of the secondary component, and it would be in fact 
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Figure 15. The relation between the orbital period-superhump 
period excess and the orbital period for dwarf novae stars. Small 
points represent known objects, from lOlech et al ] <|201lt ). The 
black dot corresponds to the position of IX Dra obtained from 
P or b and P s h- The black triangle indicate the position of IX Dra 
calculated from Pj i and P s h- Open symbols show positions of 
other ER UMa-type objects: triangles stand for ER UMa and 
V1159 Ori, and the circle for DI UMa. 



the most prominent for dwarf novae with very short orbital 
periods, i.e. ER UMa-type stars. Nevertheless, the interpre- 
tation of the observed additional humps in the case of IX Dra 
is still open and needs further studies. 

Let's get back to Tab. [5] When it comes to normal out- 
bursts and the quiescent state, the power spectrum analysis 
gave us only the information that there exists some periodic- 
ity in the data. The O — C analysis provided the quantitative 
result of a constant value P or b = 0.06646(2) d, which is in 
perfect agree ment with the valu e of orbital period of IX Dra 
estimated bv loiech et ãfl j2004Í ), P orb = 0.06646(6) d. 

With the obtained values of P s h and P or b we can draw 
conclusions about the evolutionary status of IX Dra. In 
Fig. [15] we show the position of the object in the standard 
Porb vs. e diagram for dwarf novae stars, where e is the su- 
perhump period excess defined as: 

AP P s h — Porb ,„s 

e= — h = P - ^ 

J orb J orb 

The period excess can be expressed in t erms of the mass 
ratio of the binary system, q = M2/M1 l)01ech et alJl20o3 . 
and references therein): 

e* °- 23g (7) 
l + 0.27g V ' 

Thus, the P or b vs. e relation provides us an excellent plane 
to examine the evolution of the stars. 

When we apply the formula[6]to the final set of periods, 
Psh and Porb, which were described above, we obtain the 
position of IX Dra in the region of so-called period bouncers 
(the dot in the Fig. I15[). In the same region there is located 
also another ER UMa-type star, DI UMa. 

However, it was tempting to determine the location of 
IX Dra for another prominent peak of the periodogram, 



which was presented in Sec. 14.1.21 namely P/ 4 . Although 
we did not confirm this value by the O — C analysis, its am- 
plitude in the power spectrum was slightly higher than the 
one of . The resulting position of IX Dra in the P or b vs. e 
diagram was surprising. In the case of P/ 4 , the object is 
located in the area where typical SU UMa stars are found 
before they reach the minimum period during their evolu- 
tion (the triangle in the Fig. I15[) . In the same region there 
are located also other two of five known ER UMa-type stars, 
ER UMa and V1159 Ori. 

Previous studies show that IX Dr a is located in th e 
period bouncers region in this diagram (|Qlech et al.ll2003 ). 
This is supported by our observations, but we cannot ex- 
clude the other possibility, related with the P/ 4 , which was 
discussed here. 



5 ER UMA STARS 

After 17 years since the distinction from typical SU UMa 
stars (|Kato fc Kuniaval [lÕÕol : iRobertson et ali Il995h . the 
ER UMa-type of dwarf novae is still very poorly understood. 
Until now there are discovered only five stars of this kind. 
Thus, studying each one of them is extremely important in 
order to draw conclusions regarding this class of particularly 
active dwarf novae. Because of the fact that the behaviour of 
ER UMa-type stars is so untypical, in many cross-sectional 
statistical studies of dwar f novae they are o ften s imply dis- 
count ed, see for instance iKato et al.l (|2003h and iPattersonl 
(|2011T ) . That is why particular attention shall be paid to 
investigate this class of stars separately. 

We performed a literature survey of ER UMa-type 
dwarf novae and, together with our new results, we present 
a new set of basic statistics. See Tab. [6] 

The first evident conclusion from Tab. [5] is the fact that 
the class of ER UMa-type stars is not homogeneous. It is pos- 
sible to divide them into two separate groups of objects with 
similar characteristics, which determine their evolutionary 
status (on the assumption that these extremely active stars 
obey the standard relation for typical dwarf novae stars, see 
Eq^andEJ). 

ER UMa and V1159 Ori are very alike and they form 
the first of these groups (group I). With the highest pe- 
riod excesses, they are located in the same region of the 
Porb vs - e diagram, unlike the rest of the objects. They seem 
to have normal secondaries and evolve towards the mini- 
mum orbital period. They are the only two members of the 
ER UMa class which show occasionally n egative superhumps 
l|Patterson et al.lll995l : iGao et ali 1 19991 ) . Durations of their 
superoutbursts {D s = 20 d) are significantly longer than of 
the other objects (D 3 — 10 — 15 d). 

On the other hand, IX Dra with the first version of 
parameters (see Tab. [6} is very similar to DI UMa, espe- 
cially when it comes to the period excess and location in 
the Porb vs. e diagram. Thus, they have the same evolution- 
ary status and we believe that they are period bouncers, 
i.e. more evolved objects with degenerate brown dwarfs as 
secondaries, evolving now toward longer periods (group II). 
The second set of parameters of IX Dra is possible, but less 
likely. In this case the object would be a much younger star 
still evolving towards the period minimum, so it would be 
more similar to ER UMa and V1159 Ori. 
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Table 6. Properties of known ER UMa-type stars. P or i, and P a h are the orbital and supcrhump periods. P BC and P c are the supercycle 
and cycle lengths, D s and A s are the duration and amplitude of supcroutbursts, and D n with A n are the same for normal outbursts. 
e is the period excess, and e_ is the period deficit in the case of a presence of negative superhumps. Reg. stands for region in the 
Porb vs - e diagram (see Fig. 11511 : the black triangle symbol indicates a position in the region of typical SU UMa stars before they reach 
the minimum period during their evolution, and the dot corresponds to a position of an object in the area of period bouncers. The last 
two rows show both possible sets of parameters of IX Dra, depending on the choice of the orbital period (see Sec. I4,3l for details). 
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The last object, RZ LMi, is the least known one. Its 
outbursts parameters are similar to those of other ER UMa 
stars, but RZ LMi is definitely the most active among 
them (P sc =19.07(4) d). This object shows permanent su- 
perhumps, both in superoutbursts and quiescence. An in- 
teresting thing is the fact that phases of its sup erhumps are 
coher ent even between different superoutbursts (jOlech et al.l 
2008). Unfortunately, there is no orbital period reported 
for this object, so its evolutionary status remains unknown. 
The superhump period of RZ LMi is very similar to that of 
DI UMa. Perhaps its P sc will evolve towards longer values 
over next years (like we find for IX Dra), and RZ LMi will 
resemble objeets of the group II. 

Despite the possible division of ali these objeets into 
groups I and II, there are a number of similarities between 
them. They ali have extremely short supercycles, and rela- 
tively low amplitudes of outbursts. Also, V1159 Ori shows 
a double-peaked structure within its oscillations during su- 
peroutburst, which pe ak near the end of the superoutburst 
(jPatterson et al.l 1 199a) . resembling IX Dra and its possible 
additional source of flux on the surface of the aceretion disk. 

In general, we u nderstand th e overall behaviour of 
ER UMa-type of stars. Osakij (119950 tried to reproduce the 
light curve of ER UMa, based on the standard thermal-tidal 
instability model. He increased the mass transfer rate by a 
factor of 10 compared to that expected in the standard cat- 
aclysmic variable evolution theory to obtain an appropriate 
result. With this model it is possible to reproduce normal 
and superoutbursts, together with superhump phenomenon 
during superoutbursts, which are defined as a beat between 
the precession of the aceretion disk and the orbital motion 
of the binary. 

However, it is difhcult to reconcile such an explanation 
with period bouncers. From the evolution theory we know 
that it is unexpected and unusual for the most evolved sys- 



tems to have such a high mass aceretion. A typical dwarf 
nova evolves from long (~ 10 h) to short (~ 70 min) orbital 
period due to angular momentum loss. This stage of evolu- 
tion corresponds to the mass transfer rate decrease from 
IO -8 to 10 -10 M©/yr. After the star reaches the minimal 
period, the secondary becomes a degenerate brown dwarf 
and the system starts to evolve again to longer orbital pe- 
riods. At this stage the mass transfer rate d rops to about 
10~ 1 2 M©/yr and the object becomes inactive (|Knigge et al.l 
l201ll ). Thus, the older the system, the less active it should 
be. IX Dra and DI UMa as period bouncers with brown 
dwarf donors do not fit this explanation. 

Nonetheless, the behaviour of ER UMa-type stars seems 
to be very complex. We need some improvement in the ex- 
isting models to explain also the complicated shape of the 
oscillations in the light curves. iGao et al.l (|l999f ) claims that 
"it is more likely that superhumps occasionally exist at es- 
sentially ali phases of the eruption cycles of ER UMa stars", 
not only during s uperou t bursts . which is in fact observed. It 
was sugg ested bv ISmakl (|2004 a,b) that some problems with 
the TTI model for superoutbursts might be solved by a hy- 
brid model which could combine the TTI model with the 
irradiation-enhanced mass-transfer effeets. This could im- 
prove modelling of the observed superhump amplitudes. Ef- 
feets of the irradiation-enhanced mass outflow are the most 
significant in systems with short orbital periods, so it is an 
important, probably crucial, issue in modelling ER UMa- 
type stars. 

Another important problem is the lack of high-quality 
spectra for these objeets. They are certainly needed to ob- 
tain precise orbital periods, but also to verify the nature 
of secondary stars in these systems. Ilshioka et al.l (l2007f ) 
showed that spectra of IX Dra exhibit no features of the sec- 
ondary star. It was impossible to distinguish between models 
of M2 dwarf from others. They suggested that it is possible 



14 M. Otulakowska-Hypka, A. Olech et al. 



"that the masses of the secondary stars in the ER UMa Sys- 
tems are similar, but their temperatures are higher than 
those of normal SU UMa stars". This would have a signif- 
icant influence on the evolutionary status of these objects, 
and would verify the application of the standard relation of 
typical dwarf novae stars (Eq.|6} for ER UMa-type of stars. 



6 CONCLUSIONS 

We report results of a new world-wide observing campaign of 
IX Draconis, together with our analysis of the power spectra 
and the O-C diagrams. We discovered that the supercycle 
length P sc has been increasing over the last twenty years, 
and it is now equal to 58.5 ± 0.5 d. The rate of the in- 
crease of the supercycle length is P — 1.8 x 10~ 3 . Normal 
outbursts appeared irregular during our observations, with 
the interval between two successive normal outbursts be- 
tween 3.1 — 4.1 d. We detected a double-peaked structure 
of superhumps during superoutburst, with the secondary 
maximum peaking near the end of the superoutburst. The 
mean superhump period observed during superoutbursts is 
P sh = 0.066982(36) d. We found two possible values of 
the orbital period: the first one, 0.06641(3) d, which is in 
agreement with previous studies and our O — C analysis 
(0.06646(2) d), and the second one, 0.06482(3) d, which is 
less likely. The evolutionary status of the object changes dra- 
matically, depending on the choice between these two values. 

To conclude, we updated available information on 
ER UMa-type stars and present a new set of their basic 
statistics. We emphasise the lack of spectroscopic observa- 
tions for these objects, which could provide many answers 
regarding this class of stars. 
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